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We apply a BUU type transport model for the interpretation of the di-electron 
invariant mass spectrum measured by the HADES collaboration for the reaction 
C(2 AGeV) + C. Our model incorporates the propagation of broad meson reso- 
nances emerging from the decay of bar yon resonances. 
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1 Introduction 

Di-electrons serve as direct probes of dense nuclear matter stages during the course of 
heavy-ion collisions. The superposition of various sources, however, requires a deconvolu- 
tion of the spectra by means of models. Of essential interest are the contributions of the 
light vector mesons p and u. The spectral functions of both mesons are expected to be 
modified in a strongly interacting environment in accordance with chiral dynamics, QCD 
sum rules etc. pQ. After the first pioneering experiments with the DLS spectrometer |2j 
now improved measurements with HADES |2J start to explore systematically the baryon- 
dense region accessible in fixed-target heavy-ion experiments at beam energies in the few 
AGeV region. 

We apply here our BRoBUU transport model to the di-electron data measured by 
HADES in the reaction C(2 AGeV) + C. Some features of the code are described in 
section 2, while section 3 is devoted to the presentation of simulation results and the 
comparison to data. 
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2 BRoBUU code 



The BRoBUU computer code for heavy-ion collisions developed by a Budapest-Rossendorf 
cooperation solves the Boltzmann-Uhling-Uhlenbeck equation in the quasi-particle limit £33 

dF dH dF dH dF 



a + 9F&-^= C ' tf = V(-o + [%,*))= + P* (1) 

for the one-body distribution function F(x., p, t) of a certain hadron species. This equation 
is applied to the motion of different hadron species, each with rest mass mo, in a momentum 
and density dependent mean field U. The scalar mean field U is chosen in such away that 
the Hamiltonian H equals H = \/ml + p 2 +U nr with a usually in a non-relativistic manner 
calculated potential U nr 

u" = a± + b(±Y + c ± MV Mt'rt (2 ) 

n \n J n J (2tt) 3 l + (p^) 2 

where the parameters A, B, C, r, A define special types of potentials, while n, no and Jn 
stand for the baryon number density, saturation density and nucleon distribution func- 
tion. The BRoBUU code propagates in the baryon sector the nucleons and 24 A and 
iV* resonances and additionally ir, rj, a, uj and p mesons. Details will be reported else- 
where. Different particles species (each described by a corresponding distribution F) are 
coupled by the collision integral C which also contains the Uhling-Uhlenbeck terms re- 
sponsible for Pauli blocking in the collision as well as particle creation and annihilation 
processes. The set of coupled Boltzmann-Uhling-Uhlenbeck equations is solved by using 
the parallel-ensemble test-particle method, where each distribution function is represented 

by F = 5' 3 )(x — Xj(t))5' 4 ' ) (p — Pi(t)). This method transforms the partial differential- 

8=1 

integro equations into a set of ordinary differential equations (looking like equations of 
motion) for a large number of test particles simulating the ensemble averaging process for 
the respective function F. 

Recently theoretical progress has been made in describing the in-medium properties 
of particles starting from the Kadanoff-Baym equations for the Green functions of the 
particles. In the medium particles have a finite life time which is described by the width V 
in the spectral function A of the particles. The spectral function being the imaginary part 
of the retarded propagator in the Kadanoff-Baym equations, —2ImG ret (x,p), is essentially 
defined by the self-energies S of the particle in the medium. For bosons it reads 

Av) = t^, (3) 

(E 2 - p 2 - m 2 - Re£^) 2 + \Y 2 

where ImS ret = |f and ReS rei depend on E, p and the local medium properties. The 
spectral function can significantly change during the heavy-ion collision process and can 
be simulated by an ensemble of test particles with different masses. The change of the 
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spectral function is now given by time variation of the test particle mass rrii [HI Ej • For 
bosons this additional equation reads 



_L „ (_ fleS « + ^ _D, (4) 

where the values of the selfenergy Y7 et are taken at the positions of the test particles 
i, and 5/5t stands for the comoving time derivation. The real part of the selfenergy is 
related to the mean field U. Eq. (j^J) is actually a corollary of three equations describing 
the propagation of energy, three-momentum and position of a test particle. This equation 
ensures that resonances are propagated towards their vacuum spectral functions at freeze- 
out jS]. This technique, allowing for a consistent propagation of broad resonances, is 
applied in the BRoBUU code for calculating the di-electron emission of uo and p mesons. 

Besides the propagation of broad resonances the spectral function also controls the 
production of mesons. Vector mesons (V) are essentially created by the decays of baryon 
resonances (R). The mass distribution reads 

— oc A(m v ) m v (m\ + m 2 N - m v ) X 1/2 (m 2 R , m%, m v ), (5) 

amy 

where X(a 2 ,b 2 ,c 2 ) = (a 2 + b 2 — c 2 ) 2 — Aa 2 b 2 denotes the triangular function. The reso- 
nances R is created in nucleon-nucleon collisions iViV — > NR and meson-nucleon collisions, 
MN —yR.R may more generically decay as R — > R'M into other channels (_R'=A(1232), 
iV(1440), iV(1520), iV(1535)). Resonance parameters for NN ^ NR and MN ^ R are 
fitted to available data for meson production in nucleon-nucleon and pion-nucleon reactions 
i. 

In our calculations we employ a simple form of the selfenergy: 

71 

ReS" et = 2m v Am v — , (6) 

n 

ImS- = m v (rr + ^pL=). (7) 

v 1 — v z 

Eq. © describes schematically a mass shift proportional to the density n of the surround- 
ing matter related to the normal matter density no- The effect of the collision broadening 
is given in Eq. (J2J) which depends on density, relative velocity v and the cross section ay of 
the vector meson in matter. This cross section is calculated via the Breit-Wigner formula 

= 4tt v 2 J R + 1 sT v , R Tf 

for forming resonances with masses m#, angular momenta Jr, partial widths Typt, total 
widths T^' with energy \fs and relative momentum q in in the entrance channel. In vacuum 
the baryon density n vanishes and the resulting spectral function.A mc is solely determined 
by the energy dependent width Vy c . 
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The di-electron production from vector meson decays V —>■ e + e~ is calculated by 
integrating the local decay probabilities along their paths in the collision. The subleading 
so-called direct channel — > p — > e + e~ is treated with the cross section 



/ \ m A An f a\ 2 / 4m 2 m\ t 



m p ,oAoc' 3 VM7 V M 2 (M 2 — m p ) 2 — (m p f 

with m p = iripfi + Am p n/no, m p = 775 MeV, m p = 761 MeV, T p = 118 MeV. Analog 
considerations apply to the channel np — > 

We also include into our simulations a bremsstrahlung contribution which is guided by 
a one-boson exchange model adjusted to pp virtual bremsstrahlung and transferred to pn 
virtual bremsstrahlung [T3]. Actually we use = jj<7±-^ f dqq 2 / gQ_R 2 (s) / i?2(s) and 

a pp = with <7_|_ = er pn!tot (s) s 2^i N , where M is again the e + e~ invariant mass, R 2 denotes 

the two-particle phase space volume, y/s stands for the c.m.s. energy in a nucleon-nucleon 
collision, s is the reduced energy squared after the dilepton emission, and cr pnjto4 (s) is the 
corresponding total cross section. 



3 Results and comparison to data 

We employ the above described code for the reaction C(2 AGeV) + C, where first data 
from HADES [lOj are at our disposal. In the present explorative study we are going 
to contrast simulations with and without medium modifications of p and u mesons to 
elucidate to which degree medium effects can become visible in the light collision system 
under consideration. In doing so we use fairly schematic medium effects condensed in a 
"mass shift" of Am^ = —50 MeV for the u meson. Such a shift is suggested by recent 
CB-TAPS data [11 . The use of QCD sum rules ^2] then can be used to translate this 
shift into a significantly larger shift for the p meson (dictated essentially by the Landau 
damping term); we use here Am p = —200 MeV. 

To illustrate the time evolution of the test particles we have shown in the lower part of 
Fig- ID the masses of four test particles randomly selected out of the ensemble representing 
the spectral function of an u meson as a function of the collision time. (These curves were 
obtained by artificially switching off the decay of the u mesons. Otherwise they would 
decay in a very short time in the high density region. Thus, particles with masses strongly 
deviating from their peak mass have little chance to radiate di-electrons off.) The upper 
part shows the respective local densities of the matter surrounding these particles. The 
tendency to acquire their vacuum mass starting from very different initial masses is clearly 
seen. For the test particle depicted by the dash-dotted curve which moves with a velocity 
of 0.6 c we have also shown the spectral function of the uj meson at the collision times of 
6 fm/c and 10 fm/c together with the vacuum function in Fig. EJ 

The obtained di-electron spectra are represented in Figs.H3andE] In Fig. Elthe vacuum 
parameters are employed while Fig. 0] shows the results obtained by including the above 
described medium modifications of p and u mesons. For comparison with the data the 
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Figure 1: Time evolution of the masses (lower panel) and the surrounding baryon densities 
(upper panel) of four randomly selected test particles. 



HADES filter has been applied accounting for the geometrical acceptance, momentum cuts 
and pair kinematics. The filter causes a reduction of the strength and a smearing of the 
invariant masses of the di-electrons. The result of this filtering is always shown on the 
right hand panel of the figures. 

In these figures we show different contribution to the di-electron rate. Important di- 
electron sources are 7r° and rj Dalitz decays which are proportional to the multiplicities of 
their parents. The TAPS collaboration has measured ^3] the 7r° and r] production cross 
sections of 707±72 mb and 25 ±4 mb for the same system and same bean energy. These 
values have to be compared to our calculations yielding 700 mb and 20 mb. While the 
value for pion production is in good agreement, there is a slight underestimation of the rj 
production. (Note that the presently employed cross sections rely on a global fit of many 
elementary reactions which is not optimized for the rj channel.) The Dalitz decays of p and 
uj mesons and nucleon resonances do not contribute essentially. The employed elementary 
cross sections for pp — > ppuj and pn —>■ pnuj agree with the ones in which are adjusted 
to near-threshold data for pp — > ppu. One recognizes comparing Fig. El and 0] that the 
medium modification due to collision broadening diminishes the yield of p di-electrons. 
Still an indication of a shoulder is expected from the calculations. Our calculations do not 
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Figure 2: Spectral functions for the u meson at densities felt by the test particle displayed 
by the dash-dotted line in Fig. Q at collision times of 6 fm/c, 10 fm/c and in the vacuum, 
respectively. 

show a considerable effect of a p mass shift, however this finding depends sensitively on 
our elementary p production rates employed in the calculations. Since the fine structure 
is not yet resolved in the data a conclusive decision cannot be made at present. 

4 Summary 

In summary we have compared results of a transport code to recent measurements of the 
di-electron invariant mass spectrum for the reaction C(2 AGeV) + C. The code incorpo- 
rates the production and propagation of broad resonances, in particular p and u mesons. 
Including even fairly schematically parameterized in-medium modifications of these vector 
mesons the emerging differences to results which employ vacuum parameters show up in 
a better agreement with the data. 
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Figure 3: Various sources of the di-electron invariant mass spectrum for C(2 AGeV) + C. 
Left panel: without filter. Right panel: with experimental filter and compared to HADES 
data JTZy . Vacuum selfenergies are used. 
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Figure 4: The same as Fig. \Q but with modified selfenergies of p andu mesons as described 
in text. 
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